Supersolid phases in the extended boson hubbard model 
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We present a comprehensive numerical study on the ground state phase diagram of the two- 
dimensional hardcore boson extended Hubbard model with nearest (Vi) and next nearest neighbor 
(V2) repulsions. In addition to the quantum solid and superfluid phases, we report the existence of 
striped supersolid and three-quarter (quarter) filled supersolid at commensurate density p = 0.75 
(0.25) due to the interplay of V\ and V2 interactions. The nature of three-quarter filled supersolid 
and the associated quantum solid will be discussed. Quantum phase transition between the two 
supersolids of different symmetries is observed and is clearly of first order. 
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Supersolid (SS) state 0, d, on which both diagonal 
and off-diagonal long-range order are broken, has been 
intensively discussed recently on various models [H, 0, 
S @, 0, Sj- This is partially due to the experimental 
advance of optical lattice that one day may investigate 
those hardcore boson models and the exotic supersolid 
phase experimentally. Furthermore, supersolid phase of 
spin models is also of great interests as these quantum 
spin systems has been suggested could be realized in real 
materials 0,0, H- 

The simplest hardcore boson model that includes only 
the nearest neighbor (nn) interaction, however, does not 
stabilize the supersolid phase on square lattice 

@,H- To 

induce the supersolid phase, one may relax the hardcore 
constraint to softcore 0] or to include the next nearest 
neighbor (nnn) interactions For the latter case, a 
striped supersolid (SSI) phase is found associated with 
the half-filled striped solid (QS1) phase, where bosons 
form stripes that break the x — y symmetry. This stripe 
structure allows the superfluid component to easily flow 
through the channels between stripes and therefore co- 
existence of both solid and superfluid ordering is pos- 
sible. Unlike striped solid, the hardcore checkerboard 
solid provides no pathway for superfluid component and 
no checkerboard SS of hardcore boson has been found so 
far, unless nnn hopping is included (loj . In this work, 
we present a comprehensive study on the phase diagram 
of hardcore boson hubbard model with nn and nnn in- 
teractions. A three-quarter filled supersolid that, like 
the checkerboard SS, preserves the x — y symmetry is 
found stabilized in a large parameter regime of V\ and 
V2. Associated to the supersolid is a three-quarter filled 
quantum solid, which share the same star-like occupation 
pattern (see inset of Fig[T]). For clarify, we call this solid 
and supersolid as the star solid (QS2) and star supersolid 
(SS2) hereafter. Interestingly, the supersolids, SSI and 
SS2, that possess different underlying symmetries com- 
pete in some parameter regimes in which, as we will show, 
first order phase transition occurs, in contrast to recent 
work on a similar model [13]. We tackle the problem 
with both quantum Monte Carlo (QMC) and variational 



Monte Carlo (VMC) methods which give consistent re- 
sult. A generic Jastrow wave function in VMC is able to 
generate qualitative features of all phases in QMC calcu- 
lations. 

We study the extended boson Hubbard model on a 2D 
square lattice with the Hamiltonian 
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where 6(6') is the boson destruction (creation) opera- 
tor and ^2 nn (^2 nnn ) sums over the (next) nearest neigh- 
boring sites. To set the energy scale of the problem, 
we let t = 1 throughout the paper. At half filling, the 
ground state can be a checkerboard solid (with wave vec- 
tor (tt, it)) for strong nn coupling V\, or a striped solid 
(with wave vector (7r,0) or (0,7r)) for strong nnn cou- 
pling Vi @| • For competing values of V\ and V2, how- 
ever, quantum frustration disfavors both solid structure 
and leads to the condensation of bosons instead, i.e. a 
superfluid ground state. Upon doping for large V2, as 
mentioned above, striped solid structure provides chan- 
nels of superflow so that extra bosons can form superfluid 
on top of the striped structure and leads to a striped su- 
persolid. Note that quantum effect eventually drives all 
bosons to participate the superflow, although the super- 
fluidity transverses to the stripes is much smaller Q • For 
dominating V\, on the other hand, addition bosons forms 
no condensate on the checkerboard solid because domain 
formation is energetically more favorable. As a result, no 
supersolid of checkerboard solid ordering is found. The 
phase diagram of half filling and result of doping close to 
half filling have been discussed in detail in reference [3|. 

Remarkably, when further increasing doping to p = 
0.75, our numerical calculations show a rich phase di- 
agram that contains a superfluid (SF) phase, a star 
solid phase which have finite structure factor S(Q)/N = 
Eijinnije^/N 2 at Q ? = (tt.tt), =M), and =(0,0, 
and supersolid phases of either star ordering or striped or- 
dering. The result obtained from QMC calculation with 
the stochastic series expansion (SSE) algorithm [TD] on 
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FIG. 1: (color online) Ground state phase diagram of Vi as 
a function of Vi for density p — 0.75. First (second) order 
phase transition is denoted by dotted (solid) lines. The inset 
shows the boson occupation profile (the star pattern) of the 
QS2 and SS2 phase. The cross shows a representative point 
of QS2 phase at p = 0.75 where the order parameters are 
plotted in Fig. [2] Lattice size of 36x36 and 28x28 are used 
with temperature /3 = 1/2L. 



square lattice is presented in FigfT] In QMC, the super- 
fluidity, given by p s — (W 2 )/A[3t is computed by mea- 
suring the winding number fluctuation. The calculation 
is done on scanning over different pi grand canonically to 
search for the right p that fix the density p = 0.75 for 
each coordinates (Vi,^) in the diagram. Due to particle 
hole symmetry of the Hamiltonian H , one must obtain 
the same phase diagram as in the FigfT] for p = 0.25. 
Hereafter we will focus on p = 0.75 but all discussion ap- 
plies to p = 0.25 as well. Let us now discuss each phases 
in more detail, starting with the star solid QS2. 

Inset of FigfT] shows the ordering of the QS2. The 
lattice contains four square sublattices with twice of the 
lattice constant. It is important to stress that QS2 is 
not a solid with three sublattices fully occupied and the 
fourth totally empty which naturally give p = 0.75. Our 
calculation shows that all sites have finite occupations as 
shown in the figure. As shown, two sublattices are identi- 
cal because of the x — y symmetry. A typical structure of 
QS2 has one of the sublattice almost fully occupied and 
the two identical sublattices have occupation n with the 
last sublattice has occupation ~ 2(1 — n). At V\ = 2.0 
and V2 = 6.0, the occupations on different sublattices are 
0.99, 0.37, 0.37 and 0.27 respectively. It is not surprising 
that quantum fluctuation and the gain in kinetic energy 
favor this ordering than the ordering with one sublattice 
being empty. One important feature of the QS2 phase is 
that, although the S(Q Q )/N is finite, it is rather small 
compare to those of striped solid QS1. In Figf3] we show 
the order parameters as a function of p with V\ = 2 and 
V2 = 6, a representative point (the cross in FigfT]) of SS2 
phase in the phase diagram, from which we obtain all 
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FIG. 2: (color online) QMC result of (a) boson density p, (b) 
superfluidity p x and p y , and (c) structure factor of wavevec- 
tors (7r,7r), (0, 7r) and (tv,0) as functions of chemical potential 
p — 2Vi — 2V 2 . Vi = 2.0, V 2 = 6.0 and lattice size is 28x28. 



different phases by varying p. 

A clear plateau appeared in the density curve signals 
the existence of a solid phase, in which the superfluid 
density p x ( y ) vanishes but S(Q Q )/N (FigfSJ) is finite and 
remains flat throughout the QS2 phase. Structure factor 
of all other wavevectors are essentially zero but S(Q )/N 
remains finite under finite size analysis. It is note that 
the striped solid at half filling has S(tt,0)/N ~ 0.2 (not 
shown). The small value of S(Q Q )/N in the QS2, there- 
fore, indicates the solid QS2 is rather soft. 

Away from the density p — 0.75 there are supersolid 
phases (SS2), as shown in Figf2] which is characterized 
by the same star solid ordering (wavevector Q ) as the 
QS2 state and appears on both increasing or reducing 
p, away from the QS2 phase. In this phase x-y symme- 
try is preserved such that S(n, 0)/N = 5(0, ir)/N and is 
about twice of S(tt,tt). All these peaks reduces simul- 
taneously away from the QS2 and vanish at the same 
critical point where SF emerges at larger p. This im- 
plies that this SS2 is a unique phase characterized by 
these wavevectors but not a mixture of striped phase. 
Remarkably, like the QS2 state, SS2 has one of the sub- 
lattices is almost fully occupied and does not participate 
on the superflow. Therefore the superfluid flows only on 
the other three sublattices. This self-pinning effect of 
one sublattice is rather rare and may worth further in- 
vestigation 12]. Transition from QS2 to SS2 is of second 
order as both p x (y) and S(Q )/N changes continuously 



across phase boundaries and no abrupt change in order 
parameters is observed. 

On the other hand, there is clearly a first order phase 
transition from SS2 to SSI as shown in Figf5] where all 
parameters exhibit a sudden jump at the p — 2V± — 2V2 ~ 
12.72. This discontinuous transition arises from the dis- 
tinct broken crystal symmetry of the two supersolids. 
Our VMC calculation, presented later, also supports the 
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FIG. 3: (color online) Ground state order parameters of a 
28x28 lattice at p — 0.75 for different V2. Inset shows the 
first order phase transition at V2 = 6.0. 
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FIG. 4: (color online) Phase diagram of p vs t for fixed Vi = 
2.0. The lines are guide to the eyes with dashed line (solid 
fine) stands for first (second) order phase transitions. The 
change of order parameters along the black dotted line at 
V2 = 6.0 is plotted in Fig. Awhile those along the red dotted 
at fixed p = 0.75 is displayed in Fig. [3] 



discontinuous phase transition. It is worthy noting that 
a recent study on the same model but with nnn hoping t' 
included also observes the QS2 and SS2 phases by Chen 
et al. [10]. Our result indicates that t' plays no signifi- 
cant role on the stabilization of both QS2 and SS2 phases 
which, indeed is a direct consequence of competition be- 
tween V\ and V2 interactions. Furthermore, contrary to 
our findings, Chen et al. observe a crossover from SSI to 
SS2. Whether it arises from t' is still unclear yet. 

Moving down from the cross along the dotted line in 
Fig. [U the width of the QS2 plateau shrinks as V% is 
reduced. When V> ~ 5.2 the QS2 phase disappears and 
the ground state changes continuously to the SS2. Or- 
der parameters as functions of V2 are shown in the Fig. 
[3] with Vi = 2 (dotted line in FigQ]). By reducing the 
nnn repulsion, the system gains kinetic energy that fa- 
vors superfluidity and soften the solid structure at the 
same time. Consequently, QS2 continuously changes to 
SS2 and to SF eventually at V% = 4.0 where all peaks 
of 5(Qo)/N vanish simultaneously. Note that there is 
a large parameter range where SS2 is stabilized at this 
commensurate density p = 0.75. 

More complicated situation arises for V\ < 1 where the 
nn repulsion is too weak to support the star-like struc- 
ture against the striped one. Fig. Q] shows the emergence 
of SSI at small V% < 1 within the SS2 regime. The phase 
transition between SSI and SS2 is again of first order 
because of the different broken translational symmetries 
(inset of Fig[3]). For vanishing Vi, neither SS2 and QS2 
stabilized so that only ground states of striped ordering 
is found for all fillings, consistent to the previous find- 
ings [3J]. In other words, the necessary condition for the 
appearance of the star-like quantum solid and supersolid 
is the competition between nn and nnn interactions. 

To further investigate the effect of finite V\ , we plot in 
Fig. d] the phase diagram of fixed Vi = 2 with varying 



p. The phase diagram is similar to the case of vanish- 
ing Vi (see ref.(3) except there are two new phases, SS2 
and QS2, emerges within the phase of SSI. Within this 
phase, increasing p (e.g. along the dashed line) such 
that p approaches 3/4, the star-like ordering becomes 
energetically more favorable than the striped ordering as 
discussed before and SS2 or QS2 is stabilized. Note that 
this happens only when V2 > 4, otherwise the SSI phase 
simply dissolves into SF phase upon increasing p. The 
emergence of SS2, and QS2 in the phase diagram reflects 
the interplay of nn and nnn interaction in the system. 

Now let us present our results by VMC in support of 
the star solid QS2 and SS2 found in QMC calculation. 
The wave function we used is the standard Jastrow wave 
function which is defined as : 

\y) = e-E*& Vt -' n * n '\9 ), (2) 

where |$o) = (^1=0 ^10) i s the non- interacting super- 
fluid wave function and N is the total number of bosons. 
In order to incorporate all kinds of phases in the same 
wave function, the pair-wise potential Ujj's are indepen- 
dently optimized by the algorithm proposed by Sorella 
To determine the phase diagram, we calculate the 
number of bosons in the zero momentum mode, the con- 
densate, Nk=o = b\ =Q bk=o, S(ir,0) and S(ir, tt) for the 
optimized wave function. In Fig. [5][a) three order param- 
eters are shown as a function of density for given Vi=4 
and V2—6- As density increases, the phase changes from 
SF to SS2 around the QS2 at p=0.25. With p > 0.28 
S(ir, 0) vanishes and the system becomes in SSI phase 
and QSl at p=0.5. This is consistent with the phase di- 
agram (Fig. [T|) obtained by QMC. In order to verify if 
the new SS2 is thermodynamically stable, we show the 
boson density as a function of the p which is calculated 
from the energy of adding add a particle to the system 
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FIG. 5: (color online) VMC results of (a) the condensate 
Nk=o (square), structure factor of wavevectors (71", vr) (circle) 
and (0, 71") (triangle) as functions of density and (b) boson den- 
sity p as a function of the chemical potential. Here Vi = 4, 
Vi = 6 and lattice size is 24 x 24. 
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FIG. 6: (color online) (a) VMC results of the condensate 
Nk=o, structure factor of wavevectors (it, it) and (0, it) as func- 
tions of V2 with Vi = 4 and the boson density n=0.25. The 
definition of the symbols are the same as Fig. 5. (b) The 
variational parameters Vi t j vs. m for V"2=2, 5.5 and 8. The 
solid lines are the fitting functions of Ae~^ rij . The lattice size 
is 24 x 24. 



fi = E(N + 1) — E(N). Two plateaux are found at den- 
sity p=0.25 and 0.5 which correspond to QS2 and QS1, 
respectively. Positive slope around the plateau at p=0.25 
manifests that the SS2 found is stable against phase sep- 
aration. Although there are discrepancies on the position 
of the phase boundaries, the Jastrow wave function alone 
captures the essential features and successfully generates 
all the observed phases in QMC. 

In Fig. El^a) we show the order parameters as a func- 
tion of V2 with Vi—4 and p—0.25. The three phases found 
in Fig. [3] for V2=2 are also observed here. The rep- 
resentative optimized variational parameters v(rij) for 
SF(V 2 =2), SS2(F 2 =5.5) and QS2(F 2 =8) are shown in 
Fig. E^b). The data (r > 2) is fitted to an exponential 
form Ae~ r /t with decay length £=2.0, 2.21 and 2.68 for 
SF, SS2 and QS2 respectively. The large value of v(rij) 
and £ in QS2 phase indicates the existence of a strong 
long range repulsion between bosons while the interac- 
tion is shorter ranged in SS2 and SF phases accordingly. 

In summary, we present numerical evidences for the ap- 
pearance of a quantum solid, a supersolid phase with star 
pattern and a striped supersolid at or around p = 0.75 or 
0.25. The competition between nn and nnn interaction is 
found to be important for observation of both QS2 and 
SS2. A detail study is given on the ground state phase 
diagrams by varying V% and V% as well as the chemi- 
cal potential p. The quantum phase transition between 
SSI and SS2 is appeared to be first order because of the 
abrupt change of translational symmetry. Our VMC cal- 
culation also support the QMC findings and the simple 
Jastrow wave function alone is adequate to generate all 
the phases consistent with the QMC calculations. 
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